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PREFACE 


This  report  develops  a  theoretical  means  for  predicting  ULF  near 
fields  at  long  ranges  and,  therefore,  the  maximum  range  of  ULF  com¬ 
munication  links.  Formulas  are  derived  using  a  five-layer,  stratified 
model  of  the  earth-air-ionosphere  system.  The  model  assumes  a  con¬ 
ductive  ionosphere  50  km  above  a  three-layer  earth,  where  a  resistive 
layer  is  sandwiched  between  two  conductive  layers.  The  results  are 
applicable  to  the  use  of  ultra-low  frequencies  for  special-purpose 
communications  beyond  100  km,  and  are  of  special  interest  to  the  U.S. 
Navy. 
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SUMMARY 


This  report  derives  formulas  for  ULF  near  fields  using  a  five- 
layer,  stratified  model  of  the  earth-air-ionosphere  system,  and  gives 
numerical  results  for  a  frequency  of  0.1  Hz  and  ranges  between  5  km 
and  5  Mm.  It  arsumes  a  semi-infinite,  finitely  conducting  ionosphere 
50  km  above  a  three-layer  earth,  comprising  a  moderately  conductive 
overburden,  a  resistive  layer,  and  a  highly  conductive  mantle.  It 
therefore  accommodates  the  two-dimensional  confinement  of  the  fields 
between  the  earth  and  the  ionosphere  as  well  as  possible  waveguide 
modes  in  the  resistive  layer.  However,  it  omits  waves  that  might  be 
guided  in  the  ionospheric  F-layer  and  ignores  lateral  variations  in 
the  earth's  conductivity. 

The  ionosphere  is  found  to  influence  the  signal  by  causing  its 
range  dependence  to  change  from  inverse  cube  to  inverse  square  beyond 
about  50  km.  That  result  increases  the  maximum  range  predicted  for 
ULF  near-field  communication  systems.  The  assumption  of  a  resistive 
layer  in  the  earth  between  the  overburden  and  the  mantle  slightly 
reduces  the  magnitude  of  the  calculated  fields.  The  layer's  primary 
effect  is  to  reduce  the  effective  ground  conductivity  seen  by  surface 
terminals  rather  than  to  support  a  strong  down-under-up  signal.  None 
of  our  calculations  predicts  that  magnetic  intensity  attenuates  as 
slowly  as  measured  for  signals  117  km  from  a  ULF  transmitter  on  the 
Olympic  Peninsula.  Analysis  of  the  strong  field  measured  at  that 
site  probably  requires  a  more  complicated  model  that  includes  an 
ionospheric  duct  or  lateral  variations  in  ground  conductivity. 
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I .  INTRODUCTION 


Research  on  the  terrestrial  propagation  of  long  electromagnetic 
waves  has  concentrated  on  low  frequencies  (LF,  30  to  300  kHz),  very 
low  frequencies  (VLF,  3  to  30  kHz),  and  extremely  low  frequencies 
(ELF,  30  to  300  Hz),  which  are  commonly  used  for  long-range  communica¬ 
tion  and  navigation.  Recently,  however,  the  U.S.  Navy  has  become 
interested  in  using  ultra-low  frequencies  (ULF,  <5  Hz  in  this  report) 
for  special-purpose  communications  beyond  100  km.  At  such  distances, 
the  ionosphere  and  the  earth's  structure  must  be  considered.  More¬ 
over,  ULF  wavelengths  are  so  great  that  all  points  on  earth  are 
within  the  near  field  of  aboveground  sources.  No  near-field  theory 
has  yet  been  developed  that  simultaneously  accounts  for  both  the 
ionosphere  and  possible  resistive  layers  underground.  This  report 
develops  a  theoretical  means  for  predicting  ULF  near  fields  at  long 
ranges  and,  therefore,  the  maximum  range  of  ULF  communication  links. 

Even  at  ELF,  transmission  distances  are  usually  longer  than  the 
inverse  wave  number,  A/2tt  (where  A  is  the  wavelength);  thus,  the  re¬ 
ceived  signal  is  dominated  by  the  radiation  field,  which  can  be  repre¬ 
sented  as  modes  that  propagate  in  the  earth-ionosphere  waveguide. 

Since  the  opposite  is  true  below  5  Hz,  ULF  near  fields  have  usually 
been  calculated  using  the  so-called  quasi-static  theory,  which  assumes 
the  earth  to  be  a  uniform,  semi-infinite  half  space  and  ignores  the 
ionosphere  [Kraichman,  1970].  Two  recent  analyses  have  presented 
double-layer  models  that  more  accurately  represent  the  ULF  propagation 
environment.  Fraser-Smith  and  Bubenik  [1980]  used  a  model  in  which 
a  sea  of  finite  depth  is  bounded  below  by  a  uniform  earth;  and 
Bannister  [1980]  developed  a  two-layer  model  of  the  earth.  Previous 
quasi-static  models  work  well  if  the  source  and  receiver  are  separ¬ 
ated  by  less  than,  say,  50  km — the  nominal  height  of  the  earth- 
ionosphere  waveguide  at  ULF.  However,  an  even  more  detailed  model 
is  needed  to  predict  propagation  at  the  longer  ranges  considered  here. 

One  would  expect  the  ionosphere  to  enhance  ULF  fields  at  long 
distances  by  causing  two-dimensional  geometric  spreading,  rather  than 


the  three-dimensional  spreading  predicted  above  a  conducting  half 
space.  Moreover,  a  resistive  subsurface  layer  bounded  by  two  con¬ 
ducting  layers  could  enhance  ULF  fields  by  supporting  down-under-up 
waveguide  modes  that  attenuate  more  slowly  than  near  fields.  Those 
possibilities  are  supported  by  both  theory  and  experiment.  Greifinger 
and  Greifinger  [1974]  assumed  a  uniform  earth  and  ionosphere  to  calcu¬ 
late  ULF  fields  at  distances  large  relative  to  the  height  of  the  earth- 
ionosphere  waveguide.  They  found  an  inverse-square  dependence  on  range, 
Instead  of  the  inverse-cube  dependence  that  occurs  at  short  ranges. 
Bostick,  Smith,  and  Boehl  [1977]  measured  ULF  fields  at  distances  up 
to  117  km  from  a  grounded  horizontal-electric-dipole  (HED)  transmitter 
and  found  them  to  be  stronger  than  predicted  by  an  inverse-cube  de¬ 
pendence  on  range. 

This  report  derives  formulas  for  ULF  fields  using  a  five-layer, 
stratified  model  of  the  earth-air-ionosphere  system,  and  gives  numer¬ 
ical  results  for  a  frequency  of  0.1  Hz  and  ranges  between  5  km  and 
5  Mm.  It  assumes  a  semi-infinite,  finitely  conducting  ionosphere 
50  km  above  a  three-layer  earth,  comprising  a  moderately  conductive 
overburden,  a  resistive  layer,  and  a  highly  conductive  mantle.  It 
therefore  accommodates  the  two-dimensional  confinement  of  the  fields 
between  the  earth  and  the  ionosphere  as  well  as  possible  waveguide 
modes  in  the  resistive  layer.  However,  it  omits  waves  that  might 
be  guided  in  the  ionospheric  F-layer  [Greifinger  and  Greifinger,  1968] 
and  ignores  lateral  variations  in  the  earth's  conductivity. 
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II.  ULF  SIGNALS  MEASURED  AT  LONG  DISTANCES 
FROM  CONTROLLED  SOURCE 


In  1976,  a  team  of  scientists  from  Pacific-Sierra  Research  Cor¬ 
poration,  the  University  of  Texas,  and  the  Scripps  Institution  of 
Oceanography  performed  an  electromagnetic  propagation  experiment  on 
the  Olympic  Peninsula  in  northwest  Washington.  The  experiment  was 
designed  to  send  ULF  signals  from  a  land-based  transmitter  to  a  re¬ 
ceiver  on  the  seafloor  [Bostick,  Cox,  and  Field,  1978],  but  supple¬ 
mentary  data  were  also  obtained  on  land,  east  of  the  transmitter,  at 
various  distances  equal  to  or  greater  than  the  height  of  the  earth- 
ionosphere  waveguide.  To  our  knowledge,  these  are  the  only  such 
measurements  ever  made  using  a  carefully  controlled  source.  The  land- 
side  data  have  been  reduced  but  not  yet  widely  distributed  [Bostick, 
Smith,  and  Boehl,  1977]. 

This  section  summarizes  the  landside  data,  which  contain  valuable 
information  on  the  structure  of  ULF  signals  and  the  media  in  which 
they  propagate.  Although  specific  to  the  geophysical  conditions  at 
the  measurement  site,  the  data  clarify  the  effects  of  geology  on  ULF 
signals  and  can  be  used  to  determine  the  validity  of  theoretical  models. 
The  landside  data  comprise  magnetotelluric  (MT)  soundings  at  many  sites 
on  the  Olympic  Peninsula,  and  electromagnetic  field  strengths  at  dis¬ 
tances  up  to  117  km  from  a  ULF  transmitter  that  operated  at  square-wave 
periods  between  0.5  and  10  sec.  The  transmitter  was  a  grounded  HED 
of  length  L  =  1.6  km  and  current  I  =  62.4  A  at  the  fundamental 
frequency. 

Figure  1  shows  resistivity  contours  inferred  frcm  the  MT  sound¬ 
ings.  Resistivity  is  seen  to  depend  strongly  on  depth  and  range. 

Figure  2  shows  the  magnetic  field  strengths  at  a  fundamental 
frequency  of  0.1  Hz  (dots  in  error  bars)  for  four  distances  that 
exceed  the  nominal  waveguide  height  of  50  km.  It  also  shows  curves 
that  illustrate  three  range  dependences.  The  solid  curve  plots  the 
field  produced  by  a  grounded  HED  on  a  uniform,  plane  earth.  Calcu¬ 
lated  with  the  standard  quasi-static  formula  [Kraichman,  1970]  for 


100 


Range  (km) 


Fig.  2--Measured  magnetic  intensity  versus  range 


an  assumed  effective  ground  conductivity  0 of  2.8  *  10  mho/m,  the 
solid  curve  agrees  with  the  measurements  at  ranges  up  to  60  km,  but 
predicts  too  small  a  field  at  greater  ranges.  The  two  other  curves, 
included  for  comparison,  illustrate  assumed  range  dependences.  The 
dashed  curve  assumes  an  inverse-square  falloff  with  distance  beyond 
about  60  km,  where  the  ionosphere  should  cause  a  transition  from 
three-dimensional  to  two-dimensional  near-field  spreading.  The  dot- 
and-dash  curve  assumes  an  inverse-square-root  dependence  beyond  about 
60  km,  which  would  be  expected  if  the  wave  were  guided  in  the  litho¬ 
sphere  or  the  ionosphere. 
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III.  MATHEMATICAL  MODEL 


This  section  summarizes  our  method  for  calculating  ULF  fields. 

The  complicated  equations  are  derived  in  the  Appendix,  which  also 
outlines  our  numerical  algorithm. 

Figure  3  shows  our  propagation  model,  which  represents  the  earth- 
air-ionosphere  system  as  five  layers.  Our  model  reflects  many  salient 
features  of  propagation,  including  confinement  of  fields  between  the 
earth  and  the  ionosphere,  and  guidance  of  waves  in  a  resistive  litho¬ 
spheric  layer  bounded  by  an  overburden  and  a  mantle,  both  conductive. 

By  adjusting  the  widths  and  conductivities  of  the  layers,  we  can  de¬ 
termine  the  sensitivity  of  the  fields  to  the  various  strata. 

Since  our  model  assumes  a  uniform  ionosphere,  it  must  neglect 
Alfven  waves  that  might  be  ducted  in  the  F-layer  [Greifinger  and 
Greifinger,  1968].  However,  inclusion  of  such  an  ionospheric  layer 
would  be  straightforward.  The  model  also  omits  the  effects  of  lateral 
variations  in  conductivity.  The  importance  of  such  variations  depends 
strongly  on  the  location  of  the  propagation  path.  Figure  1  shows  sub¬ 
stantial  lateral  variation  on  the  Olympic  Peninsula,  so  our  propagation 
model  might  be  better  suited  to  locations  other  than  the  site  of  the 
measurements  described  in  Sec.  II. 

The  most  common  ULF  transmitting  antenna  is  a  grounded  horizontal 

wire.  Such  an  antenna  can  be  represented  as  a  vertical  loop  comprising 

the  conductor  current  in  the  wire  and  the  return  current  in  the  earth. 

Thus,  we  assume  a  horizontal  magnetic  dipole  (HMD)  of  moment  M  oriented 

along  the  x-axis  and  centered  at  a  depth  d,  which  is  half  the  effective 

depth  of  the  return  current.  The  depth  d  depends  on  the  dipole 

length  L  and  the  effective  skin  depth  in  the  earth  6  ,  but  is  never 
3/2  8 

greater  than  6/2 
8 

The  general  equations  in  the  Appendix  are  too  complicated  to  solve 
analytically.  However,  we  obtain  useful  approximate  solutions  by  further 
assuming  the  ionosphere  to  be  perfectly  conducting  and  the  earth  to  be 
uniform.  Specifically,  the  Appendix  shows  that,  if  conductivity 


=  00 


(jj  =  oo  and  distance 


the  following  simple  expressions  give 


the  radial  magnetic  intensity  H  at  distances  much  shorter  or  greater 
than  the  ionospheric  height  h: 


H  « — — -r  cos  <j>  At/m  ,  r  <<  h,  r  <<  6..  , 

r  2irr 


(1) 


H 

r 


cos  <)> 


At/m  , 


6^  «  r  <<  h  , 


(2) 


H  « - x  cos  4>  At/m  ,  r  »  h  »  6  . 

r  2irhr  1 


(3) 


Using  M  =  lL6^//2  makes  Eqs.  (1)  and  (2)  equivalent  to  the  stan¬ 
dard  quasi-static  formulas  for  an  HED  above  a  conducting  half  space. 
That  equivalence  is  expected,  because  the  effect  of  the  ionosphere 

is  small  for  ground-level  fields  when  r  <<  h.  On  the  other  hand, 

2 

Eq.  (3)  shows  that  the  near  field  spreads  two-dimensionally  (1/r  ) 

3 

rather  than  three-dimensionally  (1/r  )  when  r  »  h.  Again,  that  re¬ 
sult  is  expected,  because  at  great  distances  the  source  plus  its 
images  in  the  ionosphere  and  earth  behave  as  an  infinitely  long  line 
of  sources. 

Equations  (2)  and  (3)  can  be  combined  to  obtain  the  following 
heuristic  formula,  which  is  accurate  if  r  «  h  or  r  »  h,  and  works 
moderately  well  if  r  w  h: 

M  cos  <J>  /r  +  2h\  _  P 

Hr  «  )  At/m  ,  r  »  6X  .  (4) 


The  coordinate  system  and  many  of  the  geometric  symbols  are 
defined  in  Fig.  3. 
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IV.  NUMERICAL  RESULTS 


This  section  presents  graphs  of  the  radial  magnetic  intensity  H 


calculated  numerically  from  equations  given  in  the  Appendix.  We  can 

easily  compute  the  other  field  components  from  equations  also  given 

in  the  Appendix,  but,  for  brevity,  omit  them. 

2 

All  calculations  assume  M  =  4n  A-m  and  cos  <|>  =  1;  distances  are 
expressed  in  kilometers.  The  results  may  be  scaled  to  any  situation 
by  multiplying  them  by  the  factor 


M  cos 
4tt 


where  M  is  the  actual  HMD  moment  expressed  in  units  of  ampere  meters 

squared.  To  compare  the  graphs  in  this  section  with  the  measured 

9  2 

fields  shown  in  Fig.  2,  use  M  =  2  x  10  A-m  in  the  conversion  factor 
(5) .  The  effective  HMD  moment  created  by  a  grounded  HED  varies  as 
l//fag,  where  f  is  the  frequency.  Therefore,  holding  M  constant  at 
4tt  A-m2  does  not  necessarily  correspond  to  holding  IL  constant. 


PERFECTLY  CONDUCTING  EARTH  AND  IONOSPHERE 


Figure  4  shows  H^  as  a  function  of  distance  for  the  idealized 
case  where  the  earth  and  ionosphere  are  assumed  to  be  perfectly  con¬ 
ducting,  semi-infinite  half  spaces  separated  by  50  km.  For  compari¬ 
son,  the  figure  also  shows  as  calculated  from  the  heuristic  Eq.  (4) 
The  transition  from  an  inverse-cube  to  an  inverse-square  range  de¬ 
pendence  at  r  «  2h  is  apparent.  The  simple  formula  given  by  Eq.  (4) 
is  accurate  if  r  <<  2h  or  r  >>  2h,  but  slightly  overestimates  the 
field  at  r  w  2h. 


DEPENDENCE  ON  GROUND  CONDUCTIVITY 

To  determine  the  effect  of  ground  conductivity  on  the  field,  we 

calculate  for  a  three-layer  model.  Formally,  we  take  D^  =  °°  and 

a  =  a  (see  Fig.  3),  and  vary  a  parametrically  while  holding  the 
—  8  8 


effective  ionospheric  conductivity  a  constant.  Figures  5  through  8 

1  _2  -3  -4 

show  a  series  of  such  calculations  made  for  a  =  10  ,10  ,10  , 

-5  1 

and  10  mho/m,  respectively.  It  is  difficult  to  assign  a  precise 

value  to  a  because  ULF  waves  interact  with  the  ionosphere  over  a 

broad  range  of  altitudes.  The  four  values  used  in  our  calculations 

-3 

span  reasonable  averages  over  a  broad  altitude  range,  with  cr^  10 

to  10  ^  mho/m  perhaps  representing  the  best  estimates  for  daytime. 

Figure  5  shows  that,  at  long  ranges,  the  field  depends  strongly  on 

effective  ground  conductivity  a  if  the  effective  conductivity  of  the 

-2  ® 

ionosphere  is  high  (a  =  10  mho/m).  If  r  exceeds  100  km,  H  is  an 

order  of  magnitude  greater  for  a  =0.1  mho/m  than  for  a  =  10  mho/m. 

8  8 

However,  as  shown  by  Figs.  6  through  8,  the  effect  of  ground  conduc¬ 
tivity  on  H  diminishes  as  a  becomes  smaller.  Figure  8  shows  that 

r  I  _5 

H  is  nearly  independent  of  ground  conductivity  if  a  =  10  mho/m. 

r  1  -1 
Figures  9  through  12  show  the  dependence  of  H  on  aT  for  a  =  10  , 

-2-3-4  rig 

10  ,10  ,  and  10  mho/m,  respectively.  The  long-range  signal  de¬ 

pends  strongly  on  0^  if  the  effective  ground  conductivity  is  high 
(Fig.  9),  and  weakly  on  o ^  if  the  effective  ground  conductivity  is 
low  (Fig.  12). 

As  expected.  Figs.  5  through  12  show  that  the  signal  is  degraded 
by  imperfect  ground  or  ionospheric  conductivity.  In  all  cases,  however, 
the  long-range  signal  is  stronger  than  predicted  by  the  standard  quasi¬ 
static  formula,  which  neglects  the  ionosphere  by  assuming  cr^  =  0. 

To  evaluate  the  possible  effect  of  a  lithospheric  waveguide  on 
the  ground-level  signal,  we  have  calculated  for  the  five-layer  model, 
using  several  assumed  values  for  the  conductivity  °f  the  litho¬ 
spheric  layer.  Such  a  resistive  layer  could  affect  the  signal  strength 
in  two  competing  ways.  First,  it  could  enhance  the  long-range  signal 
by  supporting  waveguide  modes  that  decrease  more  slowly  with  distance 
than  does  the  aboveground  near  field.  Second,  it  could  degrade  the 
signal  by  reducing  the  effective  ground  conductivity  seen  by  surface 
terminals . 

Figure  13  plots  H  calculated  for  a  5  km  resistive  layer  between 
r  -2 

a  1  km  overburden  of  conductivity  10  mho/m  and  a  mantle  of  conduc¬ 
tivity  1  mho/m — a  high  value.  This  model,  which  also  includes  an 


Range,  r  (km) 


Fig.  5--Ca1culated  magnetic  intensity  versus  range  and  effective  ground 
conductivity  a  for  aT  =  1 0~ 2  mho/m  and  f  =  0.1  Hz 
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Fig.  7 


Range,  r  (km) 


--Calculated  magnetic  intensity  versus  range  and  effective  around 
conductivity  o0  for  a,  =  1<H  mho/m  and  f  =  0.1  5"  3 
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Fig.  9--Calculated  magnetic  intensity  versus  range  and  effective  ionospheri 
conductivity  aT  for  a  =  ICH  mho/m  and  f  =  0.1  Hz 
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Fig.  13--Calculated  magnetic  intensity  versus  range  for  various  values  of 
lithospheric  conductivity  03,  a  layer  width  of  5  km, 

and  f  =  0.1  Hz 
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ionosphere  of  effective  conductivity  0^  =  10  ^  mho/m,  is  sketched  in 
the  insert  to  the  figure.  The  results  show  that  is  insensitive  to 
the  conductivity  of  the  resistive  layer.  Calculations  performed  for 
other  values  of  a^,  D^,  and  indicate  similarity  to  the  representa¬ 
tive  behavior  illustrated  in  Fig.  13. 

In  all  cases  considered,  the  resistive  layer  slightly  weakened 
the  signal  at  long  ranges,  compared  with  that  for  a  uniform  earth. 
Thus,  the  primary  consequence  of  the  layer  is  to  reduce  the  effective 
ground  conductivity  rather  than  to  enhance  the  field  through  a  down- 
under-up  propagation  mode. 


V.  CONCLUSIONS 


In  calculating  the  effect  of  a  uniform  ionosphere  and  a  three- 
layer,  stratified  earth  on  the  range  dependence  of  ground-level  ULF 
near  fields,  we  found  that  the  ionosphere  strongly  influences  the 
signal  by  causing  its  range  dependence  to  change  from  inverse  cube 
to  inverse  square  beyond  about  50  km.  Thus,  long-range  ULF  signals 
are  much  stronger  than  indicated  by  the  standard  quasi-static  formu¬ 
las,  which  apply  only  to  ranges  shorter  than  50  km  and  predict  that 
the  magnetic  intensity  decreases  as  the  inverse  cube  of  range.  That 
result  increases  the  maximum  range  predicted  for  ULF  near-field  com¬ 
munication  systems. 

The  assumption  of  a  resistive  layer  in  the  earth  between  the 
overburden  and  the  mantle  slightly  reduces  the  magnitude  of  the  cal¬ 
culated  fields.  The  layer's  primary  consequence  is  to  reduce  the 
effective  ground  conductivity  seen  by  surface  terminals  rather  than 
to  support  a  strong  down-under-up  signal. 

None  of  our  calculations  predicts  that  magnetic  intensity  atten¬ 
uates  as  slowly  as  measured  for  signals  117  km  from  a  ULF  transmitter 
on  the  Olympic  Peninsula.  Analysis  of  the  strong  measured  field 
probably  requires  a  more  complicated  model  that  includes  an  iono¬ 
spheric  duct  or  lateral  variations  in  ground  conductivity. 
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Appendix 

DERIVATION  OF  EQUATIONS 


This  appendix  derives  the  equations  that  describe  the  electro¬ 
magnetic  fields  on  the  surface  of  the  earth,  which  is  represented  in 
three  layers  by  the  model  in  Fig.  3  (p.  8). 

GENERAL  EQUATIONS 

— ► 

We  begin  by  defining  the  Hertz  potential  tt  as 

->  2->  -► 

H  =  k  ir  +  V(V  •  tt)  (A.  1) 

and 

E  =  -iwy  V  x  f  ,  (A. 2) 

■+  — ► 

where  H  and  E  are  the  magnetic  and  electric  fields,  respectively,  the 

2  2 

complex  propagation  constant  k  =  u)  -  iaivigO,  and  the  complex  re¬ 
fractive  index  n2  is 


2 

n 


(A. 3) 


We  represent  the  HED  antenna  and  its  ground  return  current  as  an 
HMD  of  moment  M  centered  at  a  depth  d  below  the  earth's  surface.  Spe¬ 
cifically,  we  write 


(V2  +  k2)TT  =  -M  , 


(A. 4) 


where 


M  =  M<$  (x)6  (y  )6  (z  +  d)ex  . 

/v 

Here,  e  is  a  unit  vector  and 
x 


(A. 5) 


M  =  (dipole  current ) (dipole  area)  . 
Following  Sommerfeld  [1949],  we  write 


■>  A  /N 

it .  =  a.e  +  p.e  , 

1  ix  1  y 


where  the  subscript  i  denotes  the  layers  (i  =  I,  0,  1,  2,  3)  illus¬ 
trated  in  Fig.  3.  It  can  be  shown,  after  considerable  rearrangement, 
that 


t*1  =  ~  f  j^A.U)  e±Z  +  B.Q)  e  ^  -  S.(A,  z)]  JQ(Xr 


)A  dA  (A. 


-cos 


oo 

/h 


(X)  e 


+  Di(X)  e 


(Xr)X  dX  , 


where  y2  =  X2  +  k2,  Rey.  >  0,  and  =  l/yi  exp  (-  yjz  +  d|)  if  the 
dipole  is  in  the  layer,  and  0  otherwise.  The  functions  A^ ,  B^,  C^, 
and  D.  can  be  found  from  the  radiation  condition  n  ■+  0  at  z  =  ±”  and 

l 

the  boundary  conditions 


ki°i(2i)  =  ki+l°ti+l(zi)  ’ 


3  3 

3z  ai  =  3z  ai+l 

z=z  ^  z=zi 


gi(zi^  =  Bi+l^zi)  ’ 


J_  B 

3z  pi 


-  A  e 

3z  pi+l 

^  8“l+l 
3x 

z=z . 

z=z . 

z=z 

Equations  (A.l)  through  (A. 13)  can  then  be  used  to  find  the  fol 
lowing  expressions  for  the  field  components  at  the  earth's  surface: 


M 


F4<X)  ‘  C0  +  D0  > 


F5(X)  -  (A0  -  B0)  -  <C0  +  D0)A  , 


F6(A)  -  A0  +  B0  . 


Application  of  the  boundary  conditions  for  the  model  shown  in 


Fig.  3  gives  the  following  expressions  for  A^,  B^,  and  D^: 


2  "Yld 

ao=^$a^ 


(A. 26) 


B0  =  Yo  V  6 


-Yld 


(A. 27) 


C0  = 


-Y0h  f  -2y  Dx  a 

A-i  e  P0^Pl  +  pP  e  +  (1  +  PjPp 

+  (1  -  P0>  e  ^  j^Aj (1  +  pp  +  A2(l  -  Pl)  e  ^  2J  jl 

r  _2yiDi  1 

*  (P^  +  PQ)  (Px  +  Pj)  e  +  (p0p-l  +  1)(1  +  p1p2M  ♦  (A-28> 

\  ~v  r  _2yidi  1 

D0  =  jA-l  e  |^pl  +  pP  e  +  p0(1  +  plpPj 

r  ~2yipi  1 

+  Aq  +  Ppe  +  (1  +  pipPj 

+  (1  -  pQ)  e  YlDl  j^U  +  pp  +  A2(l  -  Pp  e  j 

T  "2Y1D1  1 

-  (p '  +  pn)(p.  +  pl)  e  X  +  (Pnp:,  +  1)  (1  +  P,p;)  ,  (A. 29) 


e 


i  r  -YiDi  YnD,  ,  -Y^CD-d) 

L  =  - 7 -  A,  e  1  1  +  B,  e  -  —  e 

1  Y1  Y2  L  1  1  Y1 


-  A2  e 


-Y2di 


(1 


+  n)  > 


(A. 39j 


**[■ 


*2  V2  +  r3  I  *2  e 


-Y2<V®2> 


(1  +  t2)  -  A3  e 


-y3(D1+D0) 


1+D2>1 

) 


(A. AO) 


n  =  X2  e 


-2v2 


1  Yld 

A  =  —  K  e 

1  Y1 


(A.A1) 

(A. 42) 


1  _Yld 

B,  =  -7-  (1  -  K)  e 

1  Y1 


(A. 43) 


Y1  Y2D1  1  T  "Y1D1  ~Y1°1  1  -V  rd)l 

=  _1  e  2  1  ll_B  ll_J;_ell  (A  .44) 

y2  1  "  n  [  1  1  Y1  J 


Y2  Y3(Dl'hD2)n  Y2(Di+D2) 

A3  =  —  e  (1  -  T2)  e  A2  . 


(A. 45] 


2 

To  derive  Eqs.  (A. 26)  through  (A. 45),  we  assumed  »  0,  so  Yq  *=»  A. 
That  assumption  constitutes  the  quasi-static  approximation,  which  is 
valid  anywhere  on  earth  for  frequencies  below  a  few  hertz.  We  made 
this  approximation  to  save  computation  time.  The  results  presented  in 
Sec.  IV  show  the  radial  magnetic  intensity,  computed  from  Eq.  (A. 14) 
using  the  numerical  procedure  given  below. 


APPROXIMATIONS  FOR  DISTANCES  SHORT  OR  LONG  COMPARED 
WITH  HEIGHT  OF  IONOSPHERE 

_ —  -  -  -  3 

The  numerical  results  show  that  H  varies  as  1/r  if  r  «  h  and 

2  r 
1/r  if  r  >>  h.  The  transition  from  spherical  to  cylindrical  spreading 
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of  the  near  field  is  to  be  expected;  here  we  derive  analytic  expres¬ 
sions  for  that  dependence. 

As  an  example,  we  consider  the  case  of  a  uniform  earth  of  con¬ 
ductivity  o  and  a  perfectly  conducting  ionosphere.  For  that  case, 
Eq.  (A. 14)  becomes 


H  =  - 
r 


cos  4> 


A  +  tanh  (Ah) 


^jQ(Ar)A  -  ^  J^(Ar)^j  A  dA  .  (A. 46) 
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I(r) 


f  6  Ji<t>t  dt  ’ 

2r  Jo 


r  <<  h,  6^  . 


(A. 51) 


The  integral  in  Eq.  (A. 51)  is  a  standard  form  [Gradshteyn  and  Ryzhik, 
1980],  yielding  simply 


I  =  —~z  ,  (A. 52) 

2R 

2  2  2 

where  R  =  r  +  d  .  In  all  practical  cases,  r  >>  d,  so  Eq.  (A. 47) 
becomes 


H  r* — cos  <{>  ,  r  «  h,  6  ,  (A. 53) 

r  2irr 

which  is  the  well-known  form  for  an  HMD  on  a  uniform  half  space 
[Kraichman,  1970]. 

Equation  (A. 53)  is  valid  for  ranges  shorter  than  both  the  height 
of  the  ionosphere  and  the  skin  depth  in  the  earth.  It  is  most  useful 
when  the  skin  depth  is  comparable  to  or  greater  than  the  height  of  the 
ionosphere.  Another  simple  expression  applies  when  the  skin  depth  is 
small  and  6^  «  r  <<  h.  Then  we  write 


tanh  «  1 

and 


which  yield 


(A. 54) 


(A. 55) 


I  (r)  ** 


\  /  ='td/r 


Jx(t)t  dt 


(A. 56) 
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H  ^  cos  <j)  ,  <S^  «  r  «  h  . 

r  Trr 


(A. 57) 


Equation  (A. 57)  is  a  well-known  expression  for  the  field  of  an  HMD  on 
a  highly  conducting  half  space.  Note  that  both  Eqs.  (A. 53)  and  (A. 57) 
predict  an  inverse-cube  dependence  on  distance  if  the  range  is  less 
than  the  effective  height  of  the  ionosphere. 

Long-Range  Approximation 

If  r  »  h,  we  can  write 


(f)-T 


(A. 58) 


2  2  2 

(Yp  *  yV 


(A. 59) 


which  allow  Eq.  (A. 48)  to  be  rewritten  as 


00 

i_5l_  f  , 

rltVI  1 


(t)  dt  =  -  — f — r-  , 
r  1  +  y  h 
g 


r  >>  h  ,  (A. 60) 


where  we  have  used  y  d  <<  1.  Insertion  of  Eq .  (A. 60)  into  Eq.  (A. 47) 

O 

gives  the  following  approximation  for  the  magnetic  field: 


„  M  cos  <p  (  /-2i 

H  «  — ( - 

2irr  \61  +  /^2i  h 


r  >>  h  , 


(A. 61) 


which,  for  a  highly  conducting  earth,  becomes 


H  „M_cosJ 

2TThr 


r  >>  h  >>  6 


(A. 62) 
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These  equations  show  that  the  field  has  an  inverse-square  dependence 
on  range  for  distances  greater  than  the  height  of  the  ionosphere. 

Equations  (A. 57)  and  (A. 62)  can  be  combined  to  obtain  the  fol¬ 
lowing  heuristic  form,  which  works  well  at  nearly  all  ranges  for  a 
highly  conducting  earth: 


H  « 
r 


M  cos 


(^)  ■ 


«  r. 


(A. 63) 


NUMERICAL  INTEGRATION  OF  SOMMERFELD  INTEGRALS 

Equations  (A. 14)  through  (A. 19)  contain  the  Sommerfeld  integrals 


which  normally  converge  very  slowly.  When  the  pathlength  exceeds  an 
attenuation  length,  the  usual  procedure  is  to  evaluate  the  integrals 
in  terms  of  residues,  which  give  the  familiar  waveguide  modes  [Wait, 
1970],  However,  ULF  wavelengths  are  so  great  that  asymptotic  wave¬ 
guide  theory  is  inapplicable,  so  we  must  proceed  numerically. 

Bubenik  [1977]  outlines  a  numerical  method  to  integrate  Sommerfeld 
integrals  that  is  based  on  a  nonlinear  transformation  of  the  sequence 


S 

n 


c n 

I  F(X)Jm(Ap)  e"8(X)  dX  , 


n—  1,  2,  ...  , 


(A. 65) 


where  z  is  the  nth  zero  of  J  .  That  transformation,  discussed  by 
n  m 

Shanks  [1955],  can  be  quickly  computed  using  a  method  due  to  Wynn 
[1956].  We  have  developed  and  applied  a  code  that  uses  the  Shanks-Wynn 
algorithm  to  determine  S  while  keeping  the  sequence  S^,  ...,  as  short 
as  possible. 
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